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We present results on the A* (892)® production in proton-proton collisions at a beam energy of 
E = 3.5 GeV, -which is hitherto the lowest energy at which this mesonic resonance has been observed 
in nucleon- nucleon reactions. The data are interpreted within a two-channel model that includes the 
3-body production of A* (892)® associated with the A- or E-hyperon. The relative contributions 
of both channels are estimated. Besides the total cross section a{p -f p —>■ A*(892)® -f X) = 

9.5 ± 0.9lg g ± 0.7 ph, that adds a new data point to the excitation function of the A*(892)® 
production in the region of low excess energy, transverse momenta and angular spectra are extracted 
and compared with the predictions of the two-channel model. The spin characteristics of A* (892)® 
are discussed as well in terms of the spin-alignment. 

PACS numbers: 25.75.Dw, 14.40.Df 


I. INTRODUCTION 

The production of the strange vector meson A* (892) 
in proton-proton collisions was studied rather extensively 


at high collision energies m, with the lowest-energy 
data point at ^/s = 4.93 GeV [5]. No data, however. 
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exist in the region close to the production threshold 
y/sthr = 2.95 GeV. This is in contrast to the situa¬ 
tion with the kaon ground state, the production of which 
has been measured both inclusively and exclusively in a 
broad range of energies, including the very vicinity of its 
production threshold. The production of kaons and their 
excitations in proton-proton collisions is governed by the 
conservation of strangeness, so the simplest possible re¬ 
action reads 

p+p^ N + Y + K/K*i892), (1) 

where N stands for the nucleon and Y for the ground- 
state hyperons A(1116) or E(1189). The 3-body kaon 
production has been studied in detail at low excess en¬ 
ergies il HO], and it was established that, for energies 
close to the production threshold, the kaon production is 
mainly accompanied by a A(1116)-hyperon rather than a 
S(1189). It is of interest, therefore, to identify the prefer¬ 
able formation mechanism of the A'*-meson as well. 

In this paper we discuss the study of the AT* (892) + 
production in proton-proton collisions performed by the 
HADES collaboration. The deep sub-threshold K* pro¬ 
duction was analyzed in m for Ar-|-KCl reactions at 
a beam energy of 1.756 GeV. In view of future experi¬ 
ments at the Facility for Antiproton and Heavy-Ion Re¬ 
search (FAIR) exploring heavy ion collisions at energies 
of 2-8 GeV/nucleon, new data from proton-proton reac¬ 
tions are essential as reference measurements and input 
for transport models. This work complements our pre¬ 
vious studies of inclusive and exclusive strangeness pro¬ 
duction in proton-proton reactions at 3.5 GeV, namely 
dg, E(1385)+ [13], A(1405) [HUS], andpAT+A [l6|. 
The paper is organized as follows. Section H gives a brief 
information about the experimental setup. The particle 
selection and K* reconstruction procedure is described 
in Section HI. Section IV contains the obtained results, 
their interpretation within the two-channel model, and 
a discussion of the spin alignment measurement. The 
summary can be found in Section V. 

II. THE EXPERIMENT 

The experimental data stem from the High-Acceptance 
Di-Electron Spectrometer (HADES), installed at the 
SIS18 synchrotron (GSI Helmholtzzentrum, Darmstadt). 
The detector tracking system consists of a superconduct¬ 
ing magnet and four planes of Multiwire Drift Ghambers 
(MDC). The particle identification capabilities are ex¬ 
tended by the Time-of-Flight wall and a Ring Imaging 
Cherenkov (RICH) detector. The detector, as implied 
by its name, is characterized by a large acceptance both 
in the polar (from 18° to 85°) and azimuthal angles. The 
detector sub-systems are described in detail in [T7|. In 
2007 a measurement of proton-proton collisions at a ki¬ 
netic beam energy of 3.5 GeV was performed: the beam 
with an average intensity of about 1 x 10^ particles/s 
was incident on a liquid hydrogen target with an area 


density of 0.35 g/cm^ corresponding to total interaction 
probability of ~ 0.7%. In total, 1.2 x 10® events were 
collected. The first-level trigger (LVLl) required at least 
three hits in the Time-of-Flight wall in order to suppress 
elastic scattering events. 

III. DATA ANALYSIS 
A. K* reconstruction 

The Ar*(892)+-meson decays strongly, at the primary 
pp-reaction vertex, into a kaon-pion pair. The decay 
mode AT* (892)+ —>■ Ar®7r+ with a branching ratio of 2/3 
is particularly suited for the analysis, since the short¬ 
lived component of the AT®, the AT^, decays weakly into 
a 7r+7r“ pair (cr = 2.68 cm, branching ratio 69.2%). The 
considered final state of the A'*(892)+ decay is, there¬ 
fore, composed out of three charged pions, two of which 
are emitted from a secondary vertex. Therefore, as a 
first step of the analysis, we select events with (at least) 
three pions, identified by a selection on the {dE/dx )m dc- 
momentum plane (Fig. [^. The variable {dE/dx)MDC is 
the cumulative specific energy loss of a charged particle 
in the four MDG planes. 



polarity*momentum (MeV/c) 


FIG. 1. (Color online) Specific energy loss of charged parti¬ 
cles in MDC chambers as a function of the momentum times 
polarity. Solid curves show graphical cuts for the 7r+ selection. 
The Bethe-Bloch curves are shown by dashed curves. 

In the next step we consider all triplets of charged pi¬ 
ons (tt+tt+tt”) that were found in one event. Since two 
positively charged pions are available, two AT^-candidates 
are constructed for each triplet. Afterwards, to ensure 
that the A'g decayed away from the primary vertex, 
and, thus, reduce the combinatorial background, a set 
of topological cuts was applied to the tt+tt” pair that 
forms the ATg candidate. These were: i) a cut on the 
distance between the primary and the secondary vertex 
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d{Kg — V) > 28 mm, ii) a cut on the distance of closest 
approach between two pion tracks, < 13 mm, 

and hi) a cut on the distance of closest approach be¬ 
tween either of the extrapolated pion tracks and the pri¬ 
mary vertex DCA^ > 8 mm. Besides, a cut on the 
distance of closest approach between the pion track not 
associated to the if^-candidate (i.e. stemming from the 
if*-decay) and the primary vertex DCA^ < 6 mm, has 
been introduced. The application of the topological cuts 
reduces the double counting probability (i.e. identifying 
more than one itig-candidate in one event, which is kine¬ 
matically forbidden at this collision energy) to 6%. The 
invariant mass spectrum of the 7r“'"7r“ pairs that passed 
the topological cuts is shown in Fig. A prominent Kg 
signal is visible. In total, about 24 x 10^ iFg-candidates 
were reconstructed in events with at least three charged 
pions. 



M(jr+jr“) [MeV/c^] 


FIG. 2. (Color online) tt+tt” invariant mass distribution 
from events with at least three pions used for the Kg recon¬ 
struction. The legend delivers information on the fit p-value, 
the number of Kg over background, extracted mass, width, 
and signal to background ratio. 

In the next step we applied a cut on the invariant mass 
of the 7r+7r“ pairs, i.e. a 19.4 MeV/c^ interval (ilc) cen¬ 
tered at the Kg mass peak at 495.5 MeV/c^ and com¬ 
bined the pairs passing the cut with the remaining posi¬ 
tively charged pion. The resulting invariant mass distri¬ 
butions of the Kg — 7r+ system are shown in Fig. for 
the total sample (upper left panel) and separately in five 
transverse momentum bins. A clear signal of the sought- 
after K* is visible on top of the combinatorial background 
that can be divided into two classes: i) tt+tt+tt” triplets 
produced in reactions without strangeness involvement, 
and ii) non-resonant production of KgiT^ pairs. The total 
statistics of about 1700 K* allows for a one-dimensional 
differential analysis. Below we discuss the extraction of 
the transverse momentum spectra in detail; the proce¬ 
dure can be applied to any kinematical variable. The 
extraction of the raw (i.e. neither corrected for the lim¬ 
ited geometrical acceptance of the HADES detector nor 
for the efficiency of the analysis procedure) K* yields is 
carried out with fits of the experimental invariant mass 
distributions. After a careful study of the best way to 


approximate the experimental distributions, we used the 
following function 

f{M) = FpsiM)xViM;T,a) + P3{M), ( 2 ) 

where M is the invariant mass of K'g'K'^ pairs, Fps {M) 
— the factor that takes into account phase-space limita¬ 
tions, V (M;r,cr) — the Voigt function, and P3 (M) — 
a third degree polynomial that models the non-resonant 
background. The parameters of the Voigt function are: 
i) F — internal width of K* that was fixed to the PDG 
value of 50.8 MeV [TS] and ii) ct — the detector reso¬ 
lution (instrumental width) that was determined (and 
fixed as well) to be about 11 MeV by simulating a sam¬ 
ple of zero-width K*’s. The phase space distortion fac¬ 
tor Fps (M) depends on the production channel (A- or 
E-associated), so first we shall introduce the two-channel 
simulation model. 


B. Two channel model 

Due to the quite low beam energy, the three following 
3-body channels are expected to dominate K* produc¬ 
tion: 

p + p^p + A + K*(892) + , (3) 

p + p^p + i:° + K*(892)+, ( 4 ) 

p + p^ n + S+ +K*(892)+. ( 5 ) 

Here we neglect the contribution of the 4-body channels 
with an additional pion (p + p —>■ N + tt + Y -|- K*). 
They are energetically possible, but are expected to be 
suppressed in comparison with the 3-body channels. As 
it will be shown below, this assumption is confirmed 
by the experimental data. We make another simplifica¬ 
tion, namely we consider only one E-associated channel 
{pYPK*) out of two that are allowed. Up to the small 
differences in the masses of the reaction products, both 
channels have exactly the same kinematics and are in¬ 
distinguishable in the inclusive analysis of K*. Hence, 
we employ hereafter a two-channel model that includes 
the A- and the E-channels, where the latter has to be 
understood as the sum of the two isospin-splitted sub¬ 
channels. The invariant mass distributions of the sim¬ 
ulated K*'s were reconstructed for both channels of the 
model. Then, the phase-space factor Fps (M) was deter¬ 
mined, which takes into account the deviation from an 
ideal Breit-Wigner shape. As the contributions of both 
channels are not known a priori, the phase-space factor 
is constructed as the sum of the two individual contribu¬ 
tions from the A- and E-channel: 

Fps (M) = Aa X F^s (M) -b As X Fpg (M). (6) 

It was found that the fitting procedure is not sensitive to 
the exact contribution of each channel (i.e to the weights 
Aa and As). 
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FIG. 3. (Color online) Invariant mass spectra of K% — tt'*’ pairs (symbols) for the total sample (upper left) and in five 
transverse momentum bins. Solid curves are for the fits with[^ dashed curves depict the background, separately. The boxes 
deliver information on the fit p-value, the number of K* over background, mean mass and signal to background ratio. 


C. Raw spectra and corrections 


Fits of the experimental data, performed with the 
Eq. § allow to extract the raw yields (affected by the fi¬ 
nite acceptance of the detector and efficiency of the anal¬ 
ysis procedure). As an example, the raw pj-spectrum of 
the K* mesons is shown in Fig. Also shown are K* 
distributions corresponding to the A- and E-channels of 
the two-channel model. They were obtained in the fol¬ 
lowing way. A set of events (four-vectors of the reac¬ 
tions products) corresponding to both channels ^ and 
Q were simulated with the pluto Monte Carlo gener¬ 
ator |19j . A uniform population of the 3-body phase- 
space has been assumed: neither any kind of angular 
anisotropy has been implemented, nor any contribution 
from an N*- or A-resonance coupling to the K* — Y 
pair considered. Afterwards, these events served as in¬ 
put for the full-scale simulation procedure that includes 
the propagation of particles in the detector (the GEANt3 
code was employed), track reconstruction, etc. Finally, 
the simulated data sample was analyzed in the same way 
as the experimental one. All simulated curves shown in 
Fig. ID are normalized to the integral of the experimental 
distribution. Furthermore, an optimal mixture of the two 
channels has been determined by means of a x^-analysis, 
delivering the value for the relative E-channel contribu¬ 
tion of 0.4 ± 0.2, where 0.2 is the dominating systematic 
uncertainty determined by the variation of the experi¬ 
mental cuts, as will be explained below. The resulting 
mixed spectrum is shown in Fig.j^as well. A contribution 


of the 4-body channels {NttYK* final state) would pro¬ 
duce an even softer pt-spectrum as compared to the one 
generated by the E-channel and, therefore, is completely 
negligible in this analysis. Finally, we note that an analy¬ 
sis of the missing mass to the K*p final state, potentially 
more selective with respect to the A- or E-contributions, 
is not feasible due to the limited statistics. 



FIG. 4. (Color online) Raw pt-spectrum of K*’s produced 
in proton-proton reactions. Markers — experimental data 
with statistical uncertainties, curves — expectations from the 
channels (3), (4), and their mixture “0.6 x A -|- 0.4 x S”. 

As the constructed model (“0.6x A-|-0.4xE”) describes 
the data very well, we can use it to correct the raw ex¬ 
perimental yields. Due to the limited statistics of the 
experimental data we perform a one-dimensional correc¬ 
tion. The efficiency and acceptance depend, thus, on one 
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kinematical variable. We exemplify the procedure for the 
-variable; all other variables are treated analogously. 
For the purpose of the efficiency correction we pre¬ 
pare a histogram /, which corresponds to the pt-spectrum 
of the simulated K*'’s in the full solid angle, not af¬ 
fected by the limited geometrical acceptance of the de¬ 
tector and efficiency of the analysis procedure. A his¬ 
togram O corresponds to the pt-spectrum of A'*’s that 
went through full simulation and analysis chain. Finally, 
the ratio e(pt) = O/I is the efficiency histogram. Di¬ 
viding bin-wise the raw experimental pt-spectrum by the 
e(pt) histogram we obtain the acceptance- and efficiency- 
corrected spectrum. 


IV. RESULTS AND DISCUSSION 

Figure shows the acceptance and efficiency corrected 
Pt-spectrum of K*. The experimental data are normal¬ 
ized absolutely based on the analysis of the proton-proton 
elastic scattering channel in the HADES acceptance, as 
has been done for the inclusive di-electron analysis [2D]- 
To no surprise (cf. Fig.|^, the corrected spectrum is very 
well described by the two-channel model that assumes a 
40% contribution of the S-channel. 



FIG. 5. (Color online) K* transverse momentum spectrum. 
Markers — experimental data with statistical uncertainties, 
empty boxes — systematical uncertainties. The two-channel 
model (“0.6 x A-l-0.4 x E”) final state phase space distribution 
is shown by the solid curve. 


As mentioned already, a one-dimensional analysis, ex¬ 
emplified above with the transverse momentum variable, 
can be performed for any chosen kinematical variable. 
For completeness, Figs. and show corrected spectra 
for the total momentum and angular distribution in the 
pp centre-of-mass reference frame, respectively. 

Remarkably, all three selected (pt, Pc.m., and cos 0c.m.) 
projections of a three-dimensional single-particle phase 
space are well described by the two-channel model. 
No significant angular anisotropy is observed for the 
K* emission in proton-proton collisions: a Legendre- 
polynomial fit (dotted curve in Fig. does not deliver a 
better ^ as compared to the isotropic (by construction) 



FIG. 6. (Golor online) Same as Fig. [^ for the momentum 
spectrum in the pp centre-of-mass reference frame. 
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FIG. 7. (Color online) Same as Fig. [^ for the angular spec¬ 
trum in the pp-centre-of-mass reference frame. The dashed 
curve correspond to the fit with Legendre polynomials (only 
0-th and 2-nd order are used due to the symmetry arguments, 
resulting coefficients are shown in the inset). 


distribution of the two-channel model. The integration of 
the experimental pt-spectrum (Fig. allows to extract 
the total cross section of the inclusive A'*+ production: 

cr(p-bp^ A:*(892)+-bA) = 9.5±0.9;j]±0.7 pb, (7) 

where the statistical (first), systematic (second) and nor¬ 
malization (third) uncertainties are given. The system¬ 
atic uncertainty was estimated by a variation of the ex¬ 
perimental cuts (topological cuts plus selection via an 
invariant mass constraint). The values of the cuts used 
for this variations are listed in Table |T] In total, 1200 cut 
combinations have been tested. For each cut combina¬ 
tion new invariant mass fits were performed along with 
new efficiency corrections. Afterwards, from the distribu¬ 
tion of the total cross section a central interval covering 
68% of the outcomes has been identified. The borders 
of this interval—asymmetric with respect to the median 
value—define the systematic uncertainty quoted in Eq. 

The extracted total cross section value complements 
the AT* excitation function shown in Fig. in the low- 
energy region, where measurements were not available 
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TABLE I. Topological and the Kg invariant mass cut vari¬ 
ations used to estimate the systematic uncertainty. For all 
cut combinations the condition DCA^ < DCA^ has been 
demanded, reducing 1800 cut combinations to 1200. 


Observable 

Lower value 

Upper value 

Steps 

d{Kg — V) [mm] 

> 24 

> 40 

5 

d^+-^- [mm] 

< 7 

< 13 

4 

DCA^ [mm] 

> 5.6 

> 16 

5 

DCAf* [mm] 

< 3 

< 16 

6 

AM^o [MeV/c^] 

9.7(la) 

19.4(2cr) 

3 


until now. The HADES data point is consistent with the 
trend set by the measurements at higher energies. 



's- \f^[GeV] 

FIG. 8. (Color online) Energy {-/s — •^Sthr) dependence of 
the total cross section for the processes: i) pp —> A* (892) +A 
(red squares — world data [Tl|2l[4H8], red triangle — present 
work), ii) pp —> A*(892)“A (empty green crosses) [21llHS]i 
and hi) pp —>■ K^X (empty circles) f(9l ll0ll2T] and references 
therein). The solid (dashed) line is a £t to the A*(892) + 
(A+) data with f (x) = G (1 — (D/x)+)‘^, where x = 

The numerical values are C = 3.22 x 10®(1.04 x 10®), D = 
2.89(2.55) GeV, p = 1.19 x 10“^(1.16 x 10“^), u = 1.86(1.67). 

In comparison to the ground state, the K* carries spin 
one, so we proceed with the discussion of its polarisation 
properties as probed in proton-proton collisions. The 
spin configuration of the K* in the hnal state is described 
by the spin density matrix pmm' ■ The diagonal elements 
Pit) PoO) and pi_i define, respectively, the probabilities 
of the -|-1, 0 and —1 spin projections on the quantisation 
axis. The poo can be extracted from the angular distribu¬ 
tion of the decay products (A° or 7r+) in the rest frame 
of K* [22] via 

^ ^ [(1 “ Poo) + (3poo - 1) cos^(d)] . (8) 

The situation with poo 7^ 1/3 is referred to as the spin- 
alignment case, i.e. not equally probable populations of 


the ±1 and 0 spin projections. The angular distribution 
of interest is shown in Fig. Our fit with Eq. (|^ gives 
the following result: 

Poo = 0.39 ± 0.09(stat.)/:Q;Jg(syst.). (9) 

Within uncertainties our measurement is fully consistent 
with Poo = 1/3, i.e. no spin-alignment of K* is observed. 



FIG. 9. (Color online) Angular distribution of Kg in the 
A* rest frame. Markers — experimental data with statistical 
uncertainties, empty boxes — systematical uncertainties. The 
two-channel-model (“0.6 x A -|- 0.4 x E”) final state phase 
space distribution as simulated with PLUTO without any spin 
alignment is shown by the solid curve. The dashed curve — 
fit of Eq. Q to the data with poo = 0.39. 


V. SUMMARY AND CONCLUSIONS 

To summarize, we presented here the hitherto lowest- 
energy measurement of the A*(892)+ production in 
proton-proton collisions. The relative contribution of the 
channel p-|-p—i-p-l-S-l- K* (892)+ has been estimated as 
0.4 ±0.2. Within uncertainties of the experimental data, 
no deviations from a 3-body phase-space population have 
been identified in the kinematical distributions of K *, as 
well as no convincing signal of a spin-alignment has been 
observed. 

This measurement sets a baseline for future studies 
of K* production in proton-nucleus and heavy-ion colli¬ 
sions. For instance, by comparing the present data with 
previously measured by HADES proton-niobium colli¬ 
sions at the same beam energy, cold nuclear matter effects 
affecting the production of the A*’s can be extracted. 
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